† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant Nos. 11504033 and 11404290) and the General Research Fund of Hong Kong Research Council of China (Grant No. 15301014).
Molecular dynamics method is used to study the conformation behavior of a semi-flexible polymer chain confined in a cylinder channel. A novel helix-like structure is found to form during the simulation. Moreover, the detailed characteristic parameters and formation probability of these helix-like structures under moderate conditions are investigated. We find that the structure is not a perfect helix, but a bundle of elliptical turns. In addition, we conduct a statistical analysis for the chain monomer distribution along the radial direction. This research contributes to our understanding of the microscopic conformation of polymer chains in confined environments filled with a solvent.
Conformation transition of a semi-flexible polymer chain in a confined environment is ubiquitous in nature and has a variety of physical and biological applications. As a representative example, biological macromolecules such as DNA, RNA, and proteins, which are known to be confined in cells, have attracted more and more researchers’ attention.[1–7] A single DNA chain has a persistence length of approximately 50 nm, but it can be folded in a small cell with a diameter at only several micrometres. To understand the conformation of macromolecules in confined environments, a considerable number of studies have been conducted involving experiments,[8–14] theories,[15–22] and simulations.[2–7,17,23–34] For instance, Reisner et al. experimentally studied the physics and biological applications of DNA confined in nano-channels,[13] and found that these fascinating systems can be used to probe single-molecule conformation in environments with key physical length-scales ranging from 1 nm to 100 μm. Theoretically, Tree et al. suggested that the extension of DNA in a nanochannel is a Rod-to-Coil transition,[22] and showed that there exists a universal, Gauss–de Gennes regime that connects the classic Odijk and de Gennes regimes of channel-confined chains, especially for DNA in a nanochannel. Lv et al. simulated the self-assembly of double helical nanostructures inside carbon nanotubes.[25] The computational results indicated that the SWNT size and the polymer chain stiffness determine the outcome of the nanostructure while water clusters encourage the self-assembly of polyamide (i.e., PA) helical structures in a tube.
Numerous studies mainly focus on the insertion of polymer chains into hollow cylinder channels in a vacuum. The solvent is typically omitted or replaced by a simple frictional force for simplicity. We expect that the presence of solvent molecules could significantly impact this dynamic. Since the solvents have important functions in physical, chemical, and particularly biochemical processes,[35] it is important to investigate their effects on the conformation and dynamics of polymers. In particular, solvents may assist the self-assembly of new biological nanostructures.
In this paper, we investigate the confined conformation and distribution behaviors of a semi-flexible linear polymer chain confined in a cylinder channel using molecular dynamics method. In order to simulate a realistic solution environment, the cylinder channel is filled with solvent particles. The conformation behaviors under various values of stiffness of the semi-flexible polymer chains and radius of the cylinder channels are examined in detail. This work may supply a theoretical foundation for the future theories of conformation prediction and material fabrication.
Dahirel et al. have developed a new coarse-graining procedure for the dynamics of charged spherical nanoparticles in solutions,[36] which reasonably reproduces the dynamics of charged nanoparticles in suspensions. Additionally, Malevanets et al. have studied solute molecular dynamics in a mesoscale solvent.[37] They developed a hybrid molecular dynamics algorithm by combining a full molecular dynamics (MD) description of solute–solute and solute–solvent interactions with a mesoscale treatment of solvent–solvent interactions. Here, we describe the dynamics of polymer chain and solvent particles in MD method.
Our initial system comprises three parts: (i) a single linear polymer chain consisting of multiple monomers, which is the main interest of this research; (ii) solvent particles serving as a solution environment; (iii) a cylinder channel with a fixed radius formed by a regular array of particles, which represents the static confinement environment. The first two parts are confined in the third part. A periodical condition is applied in the axial direction of the cylinder (z-axis). In the simulation, reduced units are used for simplicity. The chain monomer, all of solvent particles, and the cylinder channel particles are set to with mass m = 1 and diameter σ = 1.
The bead–spring model is used to simulate the polymer chain, which comprises N monomers. The interactions between adjacent polymer monomers are described by the finitely extensible nonlinear elastic (FENE) potential
The non-bonded interactions between nonadjacent polymer monomers, solvent particles, polymer monomer and solvent particle, polymer monomer and cylinder channel particle, solvent particle and cylinder channel particle are described by the truncated and shifted L–J potential, which is used to avoid discontinuity of the potential
Moreover, the linear polymer chain stiffness is simulated by the bending potential
The number of monomer N in a polymer chain, which is also the chain length, is set to 50, 100, and 200. The bending parameter b is set to 0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, and 1000. The density ρ of mobile particles, including monomers and solvent particles, is fixed at 0.85σ−3. The radii of cylinder channel are set to 3, 4, and 5, which are all in units of σ. Due to the excluded volume interactions between mobile particles and static cylinder channel particles, the effective radii of the channel actually are R = 2, 3, and 4. To clearly describe the confinement, we define a physical parameter γ = N/(πR2), which is the ratio of chain length and the cross-sectional area. Its physical meaning is the monomer number per unit area, the larger the value of γ is, the stronger the confinement is.
Nosé–Hoover thermostat[35] and velocity-Verlet algorithm[38] are used. The time unit is τ = σ(m/ɛ)1/2, and the time step is Δt = 0.001τ. Each simulation is repeated 5 times to obtain statistically averaged results. All the simulations are conducted with the LAMMPS MD package.[38]
Under the confined environments formed by the cylinder channels and the solvent particles, the conformation behaviors of the polymer chains change considerably. Here we choose two samples to demonstrate the changes. Figure
In order to quantitatively describe the helical periodicity, a spatial correlation function[39] is studied
G(s) can be fitted with the empirical equation[39]
For a polymer chain with chain length N = 100 confined in a cylinder with R = 3, the values of G(s) are shown in Fig.
Figure
Figures
In order to understand the number of turns of the chain intuitively, we plot the average bond angle 〈θ〉 along a chain in Fig.
Figure
To determine the helical conformations in more details, the combination of cos(ϕi) and cos(φi) is used, as shown in Fig.
For radius R = 4, figure
In Fig.
We try to explain the reason of the helix-like formation by the entropy driving. Figures
We also study the distribution of the polymer chain monomers along the radial direction of the cylinder channel which are shown in Fig.
Figure
The above illustrations are all simulated under the same solution with density ρ = 0.85. Actually, we also conduct statistics about other densities. Figure
We have studied the conformation behavior of a semi-flexible polymer chain confined in a cylinder channel, filled with solvent. A novel helix-like structure is found during the simulation which is not a perfect helix structure, but a bundle of turns stuck together. The cross section of the helix structure is not circle, but ellipsoid, whose normal direction is not along the axis direction, but has an included angle from the axis. Moreover, the detailed characteristic parameters such as spatial correlation function and bond angle are calculated to verify the results. In addition, we conduct a statistical analysis for the chain monomer distribution along the radial direction, and find most of the monomers cling to the cylinder channel, except some overlapped situations. The statistical results of the helix formation percentages under various solution density and chain stiffness show that, at moderate solution density and moderate chain stiffness, there is a peak value to form the helix structure. Our work may strengthen the understanding of the microscopic conformations of confined polymer chains under solvent environment.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 | |
20 | |
21 | |
22 | |
23 | |
24 | |
25 | |
26 | |
27 | |
28 | |
29 | |
30 | |
31 | |
32 | |
33 | |
34 | |
35 | |
36 | |
37 | |
38 | |
39 | |
40 | |
41 | |
42 | |
43 |